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^ ■ We compute the topological entropy of the toric code models in arbitrary dimension at 

J> ! finite temperature. We find that the critical temperatures for the existence of full quantum 
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^ I (classical) topological entropy correspond to the confinement-deconfinement transitions in 
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the corresponding Z2 gauge theories. This implies that the thermal stability of topological 
entropy corresponds to the stability of quantum (classical) memory. The implications for 



P4 the understanding of ergodicity breaking in topological phases are discussed. 
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Topologically ordered (TO) states are ground states of certain quantum many-body systems 
that exhibit an order which does not rely on symmetry-breaking mechanism, and thus cannot be 
characterized by a non vanishing local order parameter [[l|]. They possess a ground state degener- 
acy which depends on the topology of the underlying space and which can not be lifted by local 
V~j ' perturbations of the Hamiltonian, and a pattern of long-range entanglement. 

X: 

^ . For pure states, we say they possess topological order if they span a degenerate ground space 

with a gap, and such that distinct ground states are locally indistinguishable, i.e. the reduced den- 
sity matrix over aiiy topologically trivial region does not depend on the choice of the state within 
the ground space [|2|-|4|]. This property implies a topological robustness under local perturbations, 
which has made this kind of order interesting for quantum computation [|5i 160 . This robustness 
property means that topological order is a property of a whole phase, and one is interested in some 
quantity that can label and detect this kind of order. It turns out that TO states are characterized by 
specific entanglement properties. First, they posses an area law with a finite universal correction 
[ItI-IIOII. This correction, called Topological Entropy (TE), has been shown to be constant within 
the whole TO phase of the toric code [11], and therefore can be used as an order parameter to label 
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the TO phases. It has been shown that this quantity is also characteristic of other TO phases like 
the Kitaev honeycomb model lll2l ll3ll. topological dimer phases, fractional quantum Hall liquids 
[|l4l ll50. and TO phases with finite correlation length [16], or other quantum states that are defined 
in group theoretic terms [17]. Also, these properties show up in the entanglement spectrum. To 
what extent which TO phases can be classified using the full entanglement spectrum is still an 
open problem |[l8-20]. Moreover, quantum phases break in domains of quasi adiabatic continuity, 
namely the set of those states that can be connected by means of evolution with a local Hamil- 
tonian without closing a gap [21]. If we ignore symmetry, all the non TO states are in the same 
phase as the completely factorized state. In this sense, non TO states have trivial entanglement, 
while TO states belong to different classes of non trivial, long-range entanglement [|22||. 

However, a real physical system is always found at some finite temperature T by coupling 
to a thermal environment and so does not occupy solely the ground state. If topological order 
is to be a physical phenomenon, it must therefore exist in thermal states at nonzero temperature 
too. We should then seek generalizations of the above characterizations of topological order to 
T > 0. The given definition, in terms of mass gap and ground-state degeneracy, fails to generalize 
straightforwardly as we can not speak about locally indistinguishable distinct states at T > 
simply because the thermal state is unique. Nevertheless, we can generalize this definition through 
its physical implications. At zero temperature, it implies that a TO system can support long-lived 
quantum memory. Indeed, if we add a local perturbation to a Hamiltonian of such system of linear 
size L, we must go to 0(L)-th order in perturbation theory to connect orthogonal ground states. 
Hence, the tunnelling amplitude between distinct ground states is 0(exp(— L)), in which case we 
obtain a quantum memory register with lifetime r = 0(exp(L)), which we call stable quantum 
memory. Similarly, a stable classical memory is a system in which we can reliably encode classical 
information for exponentially long times. We can generalize this viewpoint to finite temperature by 
defining a topologically ordered system to be one which supports quantum memory with lifetime 
that scales exponentially with the size of the system. As an example, the toric code in 2D [|5I] 
does not support any kind of memory, quantum or classical at any finite temperature T > 0, while 
the Ising model in 2D and the toric code in 3D both have a critical temperature below which 



classical memory is stable [|23 



m. 



The notion of TE generalizes to T > immediately [I25I1 . Scaling of the von Neumann entropy 
with subsystem proportions is more complicated then the area law at zero temperature, but one 
may take a suitable linear combination of von Neumann entropies of different subsystems and 



isolate the universal constant piece coming from structured entanglement [|9L [loll. ^ For the toric 
code in 2D and 3D, the calculation of TE at finite temperature was performed in [26] and [25]. In 
2D there is no TE at any finite temperature, just like there is no stable information. In 3D, there 
is a stable TE of completely classical origin, just like there is a stable classical memory. This 
fact strongly pushes the question whether the stability of TE and memory at finite temperature are 
always related, and if yes, why. 

Moreover, the characterization of TO as states with non trivial entanglement (NT) has recently 
been generalized to finite temperature by Hastings ^l], by considering equilibrium states that 
cannot be connected by means of a quantum circuit of finite range to a mixed state which is made 
of product states in the energy eigenbasis. Non TO states at finite temperature do possess trivial 
entanglement (FAC). 

In this paper, we investigate the important question whether these three characterizations of 
topological order are still equivalent at T > 0. We study the simplest model with TO, the Toric 
Code -in arbitrary D spatial dimensions- introduced by Kitaev ||5|], which in the low energy sector 
realizes the Z2 lattice gauge theory. We analyze the stability of quantum memory, calculate the 
TE and compare these results also with the presence of NT or FAC as indicators of topological 
order. The calculation of TE at arbitrary temperature in the thermodynamic limit is made possible 



by decomposition into contributions from the two kinds of defects [|25ll and a mapping to the Z2 
lattice gauge theory. We find the critical temperatures for the stability of TE, corresponding to the 
confinement-deconfinement transitions of the underlying gauge theory tZm- 

We find that the value of TE contains all the information about stability of quantum and classical 
memory in these models (and about the triviality of entanglement) and hence that all define the 
same notion of topological order. So, at least in the toric code models, we see that quantum 
TE means stable quantum memory, while classical TE means stable classical memory. We are 
also able to elucidate that the same physical mechanism is responsible for destruction of quantum 
memory and quantum TE. 

The toric codes examined in this paper depend on two couplings A and yU. We find there are 
two critical temperatures Tx and T^, such that the stability of memory and TE are connected as 
follows. 

• quantum memory stable and Stop{T) = StopiO) for T G [0, min(rA, T^)), and NT 

• classical memory stable and Stop{T) = StopiO) /2 for T E (min(TA, T^), max(Tx, T^)), and 
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FAC 



no 



stable memory and Stop{T) = for T G (max(TA, T^), oo), and FAC, 



where Stop{T) is the topological entropy at temperature T. 

As usual in statistical mechanics, the existence of finite critical temperature depends on the 
dimensionality of the system, where low dimensional systems are less likely to have finite- 
temperature phase transitions. As particular cases, we recover the results of Castelnovo and Cha- 
mon n25[ \26A in two and three dimensions. This precise correspondence leads us to conjecture 
that it holds in general TO systems and hence that we may define topological order at T > as 
follows: A thermal state is TO at T > if Stop{T) = Stop{0) > 0. 

The paper is organized as follows. In SectionHIl the general toric codes are described, together 
with their dualities and connections to lattice gauge theory. Statistics of defects is then used to 
analyse the stability of quantum and classical memory in the toric codes in Section Hill Decompo- 
sition of the topological entropy, duality of the toric codes and a map to Z2 lattice gauge theory 
lead to the calculation of TE in Section |lVl Our results are discussed and compared to Hastings's 
circuit definition in Section |Vl together with suggestions for further work. 



II. GENERAL TORIC CODES 



The models 

The toric codes considered in this paper are labelled by a pair {D,k), where D is the dimension 
of the lattice and A;g{1,...,-D — 1}. The toric code labelled by (D, k) will be denoted T^^'''^ 
and is defined as follows. Let A be a D-dimensional cubic lattice of linear size L with periodic 
boundary conditions and let us refer to its elementary fc-dimensional blocks as /c-cells. Let us 
denote Pa,. (A) the set of A;-cells and = the total number of 0-cells in our lattice, so that 
|Pfc(A) I = (^) A^. To obtain T^^'^\ put a spin-1/2 degree of freedom on each /c-cell and associate 
the star operator Aa = 'S'iiaedi^i with each (k — l)-cell a, where the product runs over all A;-cells 
neigbouring a, and plaquette operator Bi, = ^jg^fe^i with each {k + l)-cell b, where the product 
runs over the /c-cells contained in b. Xi, Zi are the local Pauli spin operators. The Hamiltonian of 

j-{D,k) 

i/(^-^)(A,/x) = -A 5^ 5^ P,, (1) 

aGP;,_i(A) f)ePfc+i(A) 
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where A, /i > 0. 

Aa and Bi, overlap only if a is contained in b, but then they share precisely 2 fc-cells, so that 
[Aa, Bb] = [Aa, Aa'] = [-Bfe, -Bfe'] = 0. Hence the ground state subspace is Hg = {{ip) : Aa \il>) = 
Bb lip) = 1^) Va, b}. Denote the dual lattice by A*. The algebra Ac of operators commuting with 
the Hamiltonian is generated by products of X over closed (D — A;) -chains in A* and products of 
Z over closed /c-chains in A. The algebra At = {O : O \ip) = \ip) V \ip) G Hg} is generated 
by A and B operators and hence consists of products of X over boundary (D — A;) -chains in A* 
and products of Z over boundary /c-chains in A. Hence the algebra of operators acting on the 
ground state Ag = Ad At is generated by products of X over {D — A;) -homologies of A* and 
products of Z over A;-homologies of A. The fcth Z2 homology group of the /^-dimensional torus is 
Hk{T^ , Z2) = {'1^2)^''^ ■ It is not hard to see that there is a canonical one-to-one correspondence 
between the (D — A;) -homologies of A* and A; -homologies of A, so that if we define the logical 
operators Xq,, for a = 1, . . . , (^) to be the described products over homology representatives, 
we can choose the a labels so that 

= -Z„X,, [X,,Z^]=0 ifa^/5 (2) 

The algebra of Eq.© is just the algebra of (^) independent spin- 1/2 particles, so we find 
dim-Hg = 2(fe). 

The information stored in the ground space is topologically protected since only products of 
X and Z over surfaces with nontrivial homology act nonidentically on Tig. Considering a local 
perturbation at T = 0, we would have to go to the 0(L)-th order in perturbation theory to get 
a nonvanishing matrix element between orthogonal ground states, i.e. the tunnelling amplitude 
is exponentially small in the size of the system. The situation is very different at T > where 
stability of TO depends on the energy barrier for these defects to wind around the torus. 

Duality 

The 7"''^''^^ model is exactly dual to the T''^^'^-^) model on the dual lattice, which will prove 
very useful in the following. To obtain this result, note first that each j-cell in A intersects precisely 
one {D — J ) -cell of A*, so if a G -P; (A), let a* be the corresponding {D — j)-cell in A*. Hence the 
spins naturally live on the (D — A;) -cells of A*. We now observe that if e G -P} (A) and / G Pj+i(A), 
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then 

eedf ^ rede*, (3) 

where d denotes the boundary operator. Define U to be the unitary operator swapping globally the 
X and z computational bases. Denoting A*, the analogous A, B operators on the dual lattice, 
where c G Po-k-ii^*) and d e Po-k+ii^*), we find 

UAaU^=B:,, UBbU^ = Al,. (4) 

Consequently, if 

cGPD_fc_i(A*) dePo-k+iiA") 

is the Hamikonian of the 'J'^^'^^'^^ on the dual lattice, the duality is expressed through the equation 

\i follows that any thermal expectation values calculated in the T^^'''^ and at the same 

temperature are connected by swapping A and /i since the two density matrices are conjugate. 

The unique toric code in two dimension is the well-known T*^^'^^ model with star and plaquette 
operators, which is self-dual, k also follows that there is only one kind of a toric code in 3D, since 
7^(3,1) is (jual to T^'^'^^ and so we need to go to at least 4D to find distinct models with equal D. It 
will become clear that the most interesting 4D toric code is the self-dual T*^'^'^\ 



2,2 lattice gauge theory 



In the limit A — )■ oo, the condition 



A., 



Va G Pk-i 



(V) 



is enforced on al 



11 physical states |^) projecting T^^'^^ onto the corresponding Z2 pure lattice 
gauge theory [Isl-S^ with the Hamikonian 



6ePfc+i(A) 



(8) 



and gauge constraints ([7]). Let us denote this theory by QT^^'^\ Spin flips cause excitations, which 
take the form of boundary {D — k — 1) -chains in A*. A boundary {D — k — 1) -chain C defines the 
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{k + l)-chain C* of {k + l)-cells such that Bh = -1 for b e C* . Such {k + l)-cells will be called 
"flipped cells" and C* the "flipped chain" in the following. 

Duality of QT^^^^^ and gr^D,D-k-2) for A; G {0, . . . , D - 2} can be used to show that QT^^'^^ 



for this range of /cs is a two-phase system [|29D, where the Elitzur's theorem forbids existence of 
a local order parameter. In the low-temperature phase, defects, which are necessarily at least 1- 
dimensional for A; < D — 1, are confined, and become deconfined at a finite critical temperature. 
g'-j-{D,D-i) contains 0-dimensional defects, which become deconfined already at T = 0, and thus 
the system only has the disordered phase. 

III. QUANTUM AND CLASSICAL MEMORY IN THE TORIC CODES 

General discussion 

As noted in Section HIl the primary interest in the toric codes stems from the capacity of 
the ground-state subspace to store qubits which are stable under local perturbations at zero 



temperature 113 111 . The problem of the stability of quantum memory at finite temperature is much 
more delicate and of great importance for both theoretical implications and practical reasons. In 
the view that the resilience of quantum memory at finite temperature is a property of the phase, 
one would expect that it is necessary to have a critical temperature below which memory is stable 



II32I1 . A remarkable fact of the memory encoded in the ground space of the toric codes is that their 



thermal stability can be studied through the confinement-deconfinement transition of the Z2 lattice 



gauge theory [1251 . 13111 . In this discussion, it is important to distinguish between classical and quan- 
tum memory, which are defined as follows. A qubit is prepared in a superposition |^) = q 
and coupling with a thermal bath is switched on. Let be the time scale at which the off-diagonal 
elements of the density matrix go to zero in some basis, i.e. the time at which the quantum correla- 
tions disappear. Similarly, let be the time scale when the diagonal elements change significantly, 
or equivalently when one loses even the classical probabilities. We say that our system posesses 
quantum memory if = 0(exp(L)), where L is the size of the system, and only classical memory 
if Tg = 0(1) and t^. = 0(exp(L)). Finally, the system has no memory if both r^, = 0(1) in the 
size of the system. 

The ground state subspace of T*-^'^^ is isomorphic to the Hilbert space of (^) 2-level systems, 
and its algebra of logical operators is generated by X^, Z^iox a G {1, . . . , (^) }, which satisfy 
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(l2l). The different a sectors are equivalent and independent and we will restrict to a = 1 in the 
following and drop the indices. Let us choose a basis {|0) , |1)} for this tensor factor of "Hg such 
that 2 |0) = |0), 2\\) = - X |0) = |1) and X |1) = |0), prepare qubit in a superposition 
= Co |0) + ci |1) and turn on coupling with a thermal bath at temperature T. The coupling 
is assumed local and so T > will produce local defects of Aa and Bb. Quantum memory is 
destroyed when the thermal defects can change the eigenvalue of either X or Z, while classical 
mernory is preserved when the eigenvalue of either X or Z is robust under the thermal defects 



I133i 



3411 ■ We have seen that X is a product of X over a. {D — /c) -homology in A* and Z a product 
of Z over the dual A; -homology in A. Hence the eigenvalue X is fragile when the defects can 
wind around the torus. In turn, the deconfinement of these defects can be seen as the deconfined 
phase of the corresponding gauge theory, i.e. the deconfined phase of QT^^'''^ with coupling /i. 
Similarly, the eigenvalue of X is fragile when the Aa defects can conspire to produce surfaces 
with nontrivial homology. The X fi duality tells us this happens precisely when ^7~(^'^^'^) 
with coupling A is deconfined. We can now summarize these results in the following table 

quantum memory ^ both gr^D,D-k) ^j.^^-^ gj-{D,k) ^j^^^^ confined 
only classical memory ^ either gT^^^^-''\T/ \) or GT^^^'^^T/ fx) deconfined 
no memory ^ both gr^D,D~k) (jj^^ g-yiDM) (5-/^) deconfined 

TABLE I: Quantum and classical memory in T^^'^"* 



Examples 



The lattice gauge theory QT^'^'^\T / fi) has phase transition at T/fi = 0, and so the 2D toric 
code has neither quantum nor classical memory at any finite temperature Issl]. QT^^''^\T / \) also 
deconfines at T/A = 0, but QT^^'^\T/fi) has a nontrivial low-temperature phase, so that the 
3D toric code loses quantum memory at T = 0, but preserves classical memory up to a finite 
critical temperature, which is proportional to the coupling of the plaquette operators. The simplest 
toric code with quantum memory at finite temperature is the T^^''^^ model pllj . whose stability is 
controlled by QT^^'^^T / /j) BM- 
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IV. TOPOLOGICAL ENTROPY 



Definition 



In this section, the topological entropy of T^^'^^ will be calculated for general couplings at any 
T in the thermodynamic limit and shown to correspond precisely to the behaviour of quantum and 
classical memory as discussed in the previous section. 

At T = 0, topological entropy Stop{^) is defined as the universal part of bipartite entanglement 
entropy Se, which does not scale with the subsystem boundary 

Se = a\d\-StoM- (9) 



As opposed to a, it is robust under local perturbations [lllh and also Stop{^) > 0, so that it repre- 
sents an order of entanglement. 

At T > 0, the scaling of Se is more complicated and we need to extract the universal term 



Stop{T) by taking a linear combination [l9|, 



lOD 



Stop{T) = J2^i^)Se{T,a), (10) 

i 

where a{i) are the signs of partitions Ci U Di = Pfc(A), which are chosen so that the linear 
combination of the bulk and boundary chains of Ci of any dimensionality are zero, and we are left 



wit 



hti 



le topological contributions only. Here we generalize to arbitrary D the clever construction 
of |l25]. For a general T^^'^\ we will define the Cf\ i G {1, . . . , 4(D — 1)} by induction on D 
as follows. For D = 2,we choose the four partition as in Figure[Tl with cr(2) = cr(3) = — o"(l) = 
— cr(4) = 1, where the outer, inner squares in C^^ have sidelengths a, a/3 respectively. Having 
constructed all Ci^\ we define C|^^^'' to be the {D + 1) -dimensional cube of side a missing a 
{D + 1) -dimensional cube of side a/3 from its middle. Analogously to FiglH we define C[^^^^, 
^4D~2^ as the upper and lower two-thirds of C[^'^^\ upper and lower meant in the last dimension, 
and Ci2-3^ = CS>-1 n Ci^t'l For z e {1, . . . A{D - 1)}, construct = ^ x I,/,, 

where the interval la/s of length a/3 is added in the last dimension so that C^^'t^J = C[^'^]^. 
Finally, choose a{i) 

{-1 ifi = 0, 1 mod 4 
(11) 
+ 1 ifi = 2,3 mod 4 

As an example, for D = 3 and D = A we have the partitions of Figl2l and [H The fourth 

dimension is represented by the green segments in Fig. [3l while the blue and red are 3-boundaries 
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FIG. 1: Partitions for the i:> = 2 model 

living in the first three dimensions. This choice of partitions was motivated by the requirement 
that the signed bulk and boundary chains add up to zero. Indeed, they clearly do for D — 2, and 
hence by induction on D, 

4(D-1) 

(12) 



4(D-1) 



since C^^+^^ = C^"^' x 1^/3 for i e {1, . . . , 4(L> - 1)}. Moreover, it follows straight from the 
definition that 

4D 

y: a{{)cr'^^o, (13) 



i=4:D-3 

and hence the full set of partitions gives zero net bulk chain 



41? 



0. 



(14) 



The same argument also works for the non oriented boundaries with signs a{i). When choosing 
the partitions, we also required that their collection is symmetric under the exchange Cj ■<-)■ D^, 
besides the global torus topology. For example, C[^^ has 2 connected components and so in D 
dimensions, we are forced to introduce the i = 4{D — 1) partition, where D^^_i) also has two 
connected components. Finally, if we also demand that for each o? e {0, . . . , L> — 1}, there is an 
i such that Cj has nontrivial homology of dimension d, our choice is a very natural one. In the 
thermodynamic limit, we scale both L, a — > oo. 

Zero temperature 



Let us first calculate Stop at zero temperature as this result forms the core of the calculation 
at T > 0. Define the groups of spin flips G = {Aa\a e Pk-i)-, Gi = {g e G\gDi = Izj^} 
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FIG. 2: Partitions for the D = 3 model 





FIG. 3: Partitions z = 1, 5, 12 for the i:* = 4 model 



and Hi = {g E G\ gc, = ^d}, where go, denotes restriction of g to Di and is the identity 
transformation on Di. Since any ground state is a uniform superposition over the group G, the 
entanglement entropy of partition i is O 

IGI 



Hence the TE is 



Si^''^^ = log 



4(D-1) 



I G,; 1 1 Hi I 



(15) 



(16) 



i=l 



because there is an equal number of positive and negative partitions. Now we turn to the calculation 
of and \Hi\. The contribution from the global topology of A to log (| if, |) is the same for all 
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i and therefore cancels in (fT6l) . It will therefore be omitted in some of the following equalities, 
which we will write as = for that reason. Gi consists of those closed {D — /c) -chains in C*, which 
are also boundaries of (D — k + l)-chains in A*. The elements of Gi which are not boundaries of 
chains in Gi must arise from nontrivial (k — l)-homologies of D*, so that 

log(|G,|) = hg{\H,.,{D*)\) + \og{\BD-k{G:)\), (17) 

where Hd, Bd are the Z2 homology, boundary groups respectively. Let Xd and Zd be respectively 
the groups of all, and closed Z2 d-chains in Gi. In other words, Zd = kei{dd), Bd = im(dd+i), 
where dd : Xd^Xd-i is the boundary homomorphism. Then the Z2 homology groups are defined 
by 

Hd = Zd/Bd. (18) 
By the first isomorphism theorem applied iteratively to 5^ for d G {D — k + 1, . . . , D},we obtain 

k 

log2(|G,|) = b,^,{D*) + - |P^_(,_,)(C;)|], (19) 

i=i 

where bj is the j-th Betti number. Similar analysis holds for \Hi\ and we arrive at the result 

= - V a(z) |6,_i(A) + + 
log 2 ^ 1 

k ^ 

+ J2i-^yibD-ik-,)iC.^ + bD-ik-j){D.i)] \ , (20) 
j=i J 

where the terms containing the j-cell numbers \Pj\ subtracted since the bulks of signed partitions 
add up to zero, and the stars were dropped since Gi, Di have the same topology as G*, D*. Our 
task is reduced to calculating the Betti numbers of Cj and Di. Define reduced Betti numbers 
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b'o = bo — 1 and 6- = bi for i > 0. They are 



b[{C,,) 



1 if i = j 



otherwise 




1 ifi = D-j -1 
otherwise 

6:(A,+3) = 0. 



(21) 

(22) 

(23) 
(24) 

(25) 
(26) 



Substituting these values into (|20l) . we find that 



5L^'')(0) = 21og2, 

independently of D and A;, which also reproduces the known results for 2D and 3D [|25 
we separated the contribution from and we would have found 

4(D-1) 4{D-1) 



(27) 
.a^]- Had 

(28) 



i=l 



i=l 



Finite Temperature 

Let us now proceed by calculating sI^\t, A, jj,) for general T. To perform this calculation, 
we will use the key property found in ll25ll . where it was shown that the entanglement entropy 
in general toric codes decomposes into a sum of two terms, with the first coming from the star 
operators and depending only on T/A and the second from plaquette operators depending on T//i. 
Thus we may write 



(29) 
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where the reason for notation Qi^'-'\ Ri^'-''' will become clear shortly. TE is linear in von Neumann 
entropies, so the same factorization applies 

sIo/\t, a, /.) = Qi^f-'\T/X) + <'^Ht//.). (30) 

We have seen in Section HI] that T^^''^^ with couplings (A,;u) is equivalent to the 7"^^'^"'^^ with 
couplings (yU, A) on the dual lattice. Since Cj, Di have the same topology as C*, D*, the same 
duality holds for TE 

^f5'')(r, A,/i) = ^g'^-'^HT,/", A). (31) 
It follows that Q[^p^^ and R^op^^ can be chosen to coincide 

Qiof{x)=RSf{x) (32) 

and we can write 

Si:^p'\T, A, fi) = Q\^f-^\T/\) + Q\^J'\T/^^). (33) 

'top 



We have thus reduced our problem to finding the functions Q\op'' {x) for j G — 1}. The 

zero-temperature result tells us that we must have Q^^''^^ (0) = log 2. The A — )■ oo limit of equation 
(1331) then yields 

QS'^ {T/^^) = lim Sl!^''^ (T, A, /i) - log 2. (34) 

But limx^oo Sj:l^p''\T, X, fi) is nothing but the topological entropy in the gauge theory 
QT^^'^'^T/fi) since A — )■ oo imposes the gauge constraints Aa = 1. 



Topological entropy in the gauge theory 

In order to compute Sgfgp\T/ jj) = limA^oo S[^p^\t, X, jj), introduce the projection operator 
P onto the states of the gauge theory 



The thermal density matrix of our model is 

" (36) 



\ 6GPfc+i(A) / 
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where (in the computational z-basis) 



Z = Tr 



exp I J2 P 



\G 



\G 



(37) 



is the partition function of QT^^'^\T/^), where F is the group of all spin flips on Pk{A.) and 
^if) = (/I ZlfeePfc+iCA) I/) is the total plaquette magnetization of /. Let us also define the 
groups Fi = {f e F\ fo, = IdJ and Ei = {f e F\ fc, = IcJ that act nontrivially only on d, 
Di respectively. The reduced density matrix of subsystem Cj is then 

E.eG. e-/^^*^(/) (^1/) (/Die. 

P^iT/^I) = , (38) 



where 0]^ denotes the projection of operator O onto the Hilbert space of Cj. Let us use the 
replica trie to find the entanglement entropy of Ci in the gauge theory 



d 

dn 



(39) 



n=l 



The trace of the n-th power of p" is found from (|38l) to be 



Tr [pUT/p 



(40) 



fU-JnGF 



where 



fn 



Ci 



(41) 



1 iff /lie. = 
otherwise. 

We can trivialize the delta constraint through substitution /^ = /e^, where / G Fj and G E, 

\n In— 1 



Tr [pr(r//x)] 



\G, 



n— 1 



Yl Yl 

f&Fi ei,...,e„&Ei 



-/3/iEM(/e, 



m=l 



where 



<l^if) = E 



(42) 



(43) 
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Expression is ready to be used in the replica trie 091 ). with the result 

sT\t/i^) = - log(IG.I) + \og{\G\Z) - ^ E ^^(/) 



(44) 



Observe that defines a probability distribution 



|G|Z 



(45) 



on the group Fj. Indeed, Xl/Gi^; Qiif) = 1*^1^- Taking now the linear combination over the signed 
partitions CiU Di, the i-independent terms cancel and we arrive at 



4(Z)-1) 4(D-1) 
i=l 1=1 



4(D-1) 



log 2- J2 



i=l 



X^Pi(/)logPi(/) 



(46) 



where (l28l) was used in the second equality. The first term is just the contribution of \Gi\ to the 



zero temperature result (fT6l) . Returning back to Q[^''\t/h), we find from equation (|34l) 



4(D-1) 
i=l 



(47) 



Isolating gauge redundancy 

Expression (1471) still includes some gauge redundancy. Define the groups of spin flips Q = 
G X {Xa\a = 1, . . . , (^)), = {51 G ^1 5fD. = 1a:} and = {5- G ^| gc\ = IcJ analogous 
to G, Gi and //j, but this time also containing the noncontractible flips. Magnetization M(/) is 
invariant under f t-^ gf for g E Q, so that for f E Fi 



\G\Z 



E 



(48) 



where Ei = Ei/T-Li. Moreover, if we define the group Ki = {f E Fi\3e E Ei : fe E G},we find 
thatpj(/) only depends on the coset [/] G Fi = Fi/Ki. Note that \Ki\ = IGl/lT-Li], and hence 



5^J9.(/)l0gp.(/) = log(|?/.|) -log(l^l) + 5^P,(/)logp.(/), 

where pi{f) is the non-redundant probability distribution over Fi given by 



(49) 



-2/3M-I'{/e) 



(50) 
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(51) 



where $(/) is the number of flipped {k + l)-cells in configuration / and 

F = F/Q. Elements of F are precisely the physically distinct configurations of the gauge theory, 
and elements of Fi are those configurations of (k + l)-cells inside Ci which can be extended to 
consistent global configurations of the gauge theory, i.e. to boundary (D — k — l)-chains in A*. 
Similarly, one should regard Ei as containing distinct ways to complete the boundary {D — k — 1)- 
chains inside D*. We have now completely removed the gauge redundancy from our expressions. 
Observe that by (|28]) 



4(D-1) 4(_D-1) 

J2 cx{t)\og{\n4) = - Yl ^(01og(|i/.|)=log2, 

1=1 i=l 



(52) 



and consequently 



Q 



(D,k) 
top 



4{D-1) 

(T//i)= log 2- J2 ^(^) 



1=1 



YviWogpi^f) 



(53) 



In other words, the gauge redundancy contributes the factor log 2 to TE. This factor can be identi- 
fied with the \Hi\ part of (fT6l) . 

Phase transition 



At this point, we are finally ready to show the mechanism that produces a phase transition at 
a certain temperature T = Tc for the topological entropy. To this end, we investigate the second 
term in Eq.(l53l) and show that it leads to a phase transition in q[I^p''\t / fi) . First notice that at 
T = 0, only / = 1 G Fj produces nonzero probability pi, since all other configurations contain 
flipped (k + 1) -cells. But thenpj(l) = 1, and so indeed 

Ql5')(0)=log2, 

as required by self-consistency. Let us rewrite the square bracket in (|53l) as 



E P^(f) ^^sMf) = E ^e-^^^'^^/^) log 



/6F 



^ W 

e'eEi 



-2/3/x$(/e') 



-2/3Ai*(/e") 



(54) 



(55) 
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where we used F = Fi x Ei and the fact that mukiplication hy e ^ E Ei produces a mere 
permutation within E^. All the dependence on i is now inside the logarithm and thus 



C''^(T//i)=log2-5^r(/)log 



s-{f) 



where 



and 



r{f) 



1 



-2/3m*(/) 



(56) 



(57) 



s±{f) 



exp(-2/x/3$(/ei)) 



(58) 



n 

i\(7(i)=±l 

Let us suppose that QT^^'''\T / fi) is in the confined phase. Big membranes of defects are 
suppressed, i.e. r(/) is exponentially small in the defect extent. In particular, defects can not 
detect the topology of Cj. When we expand s±{f) for fixed / G F, we can see it is a sum over 
overlapping configurations of the gauge theory, differing from / only in DiS. For / = 1, the 
allowed multiplicity of a defect at a given b E Pk+i{A) in a term of s±{f) is equal to |{i|cr(i) = 
±1 A 6* G Po-k-iiD*)}\, for the plus and minus sign respectively. Since all contributing defects 
are local and the partitions satisfy the chain equation J2i\a(i)=i ~ '^j\a{j)=-i ^^^^ 
there is one-to-one correspondence between contributing terms in s+ (/) and s^{f), and hence that 



sM) = S-if) 



(59) 



(60) 



for any f E F. Therefore, for T/ fi < {T/fi)crit 

Qiop'\T/fi) =\og2. 

On the other hand, when QT'^^'''^ (T/ /i) is deconfined, T//i > (T/ jj,)crit, a typical configuration 
will contain many topological defects, i.e. (D — k — l)-branes (which are boundaries). By (l24l) - 
(|26|) . the only D*s with nontrivial {D — k — l)-homology are those with i = Ak, Ak + 1, where 
the later only exists iik < D — 1 (the global toric topology does not count as nontrivial here as it 
contributes to all bipartitions in the same manner). Let / G F be a membrane which wraps around 
Cl^, i.e. a nontrivial (D — k — l)-homology in Dl^. Then / G Ei precisely for i < Ak, hence 



\{i\a{i) = lAlEEi}\ = \{i\a{t) = -1 A / G Ei}\ + 1. 



(61) 



When k = D — 1, every term in after expanding the product over i\(T{i) = 1, can be 

obtained by twice as many ways as the same term in This is because Dl^jj_-^^ must contain 



19 



an even number of topological defects (i.e. an even number of point-like defects inside the smaller 
hypercube), whereas any parity is allowed for i < 4(D — 1). For general k, remember that above 
Tc there is an infinite number of defects in the thermodynamic limit, and therefore on average 
again there are twice as many ways of distributing the topological defects among the positive 
partitions with respect to the negative ones. Moreover, the variance of this distribution is zero in 
the thermodynamic limit i47ll . Hence, in the thermodynamic limit, above Tc, we have 

s-if) 

for any / G F. To illustrate this point, consider the T^'^'^\ where the defects are 0-boundaries, 
i.e. pairs of points virtually connected by a line. Terms of s_(l) containing an odd number of 
endpoints inside the small square must all come from i = 1, whereas in they can come 

from both i = 2,3. Similar reasoning applies when f ^ 1. Consequently, if T/yU > (T/ ij)criu we 
obtain 



^r{f) log 



and the two terms in (|56l) precisely cancel. Hence 

log 2 if gr^^^^'^ (T//i) confined 

iigT^^^^\T/n) deconfined, 



log 2 (63) 



C''Hw= r^"" . (64) 



Qi!^''\T/fi) experiences a phase transition at the same T/fi as the lattice gauge theory. We 
can finally substitute into (133]) to find the behaviour of TE as in Table UIl The similarity with Table 
Uis striking. 

Si^/\T, X,iJ.) = 21og2 ^ both gT^^'^-^\T/\) and gT^^'''\T/ 1^) confined 
5^5*^^ (T, A, fi) = log 2 ^ either gr^D,D-k) ^^/^^ gq-{D,k) ^^/^^ deconfined 
S\^p^^ (T, A, /i) = ^ both gr^D,D~k) ^^/^^ gj{D,k) ^^/^^ deconfined 

TABLE II: Topological entropy in T^^'^) 



V. DISCUSSION 

We have shown that the topological entropy is a good order parameter for topological order in 
general toric codes at finite temperature. If we denote t^nf'' = {T/n)crit the critical coupling in 
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the QT^^'^^^T/ fi), we can conclude that the topological entropy of 7"^^'^-' experiences two phase 
transitions at Ta = A tl,^if~^^ and = /x t^^if \ such that 



2 log 2 forT < min(TA,T^) 

log 2 for min(TA, T^) < T < max{Tx, T^) (65) 
for max(TA,T^) < T. 



^crit = Oif/c = D — 1 and t).^-^ ' > otherwise. Thus in models with k = 1,D — 1, TE 



t^f ) = if A: = D - 1 and t^Sf 
has its maximal value So = 2 log 2 only at zero temperature. The well known T"*-^'^-* model is 
the only one for which Stop{T, X, fx) = for all T > 0. Models with 1 < k < D - 1 have 
a nontrivial low-temperature phase, in the sense that SJJ''^(T, A,/i) = So for all < T < T„ 
Tc = min(TA, T^) > 0. Such models only exist in D > A, the simplest example being the Kitaev's 



four-dimensional toric code 1311] (T'-'^'^-'). 

The temperature dependence of Stop follows the same pattern as the properties of quantum and 
classical memory in our systems. Indeed, if we combine Tables U and UIl we discover the promised 
connection between TE and robustness of memory in the toric codes 

quantum memory <^ s[^p^\T, A, /x) = 2 log 2 
classical memory <^ s[^p^\t, A, fi) = log 2 
no memory <^ s[^p^^ {T, \, fi) = 

TABLE III: Memory and Stop in T^^'^^ 
In the introduction, we mentioned the recent formulation by Hastings for topological order at 



T > [|27l] . A thermal state peq has TO if it can not be transformed arbitrarily close to a 'classical' 
state by means of local unitaries, even if we are allowed to tensor in additional local degrees of 
freedom. By a classical state is meant a thermal state of a local Hamiltonian which is diagonal in a 



product basis. In [27], it is argued that a thermal state is topologically ordered if one can efficiently 



perform quantum error correction, in the sense of PJQ, since one can then thicken the logical 
operators X, Z while preserving their algebraic properties. These then prevent local unitaries 
from conjugating the state peq arbitrarily close to a classical state. In T^^'^\ efficient correction 
of defects of both kinds is possible if and only if both gr^D,D-k) j^^^^^, gq-{D,k) i .^^ 
the confined phase. Hastings's circuit definition, the presence of stable quantum memory and the 
behavior of TE all agree in predicting topological order in the toric codes: 
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A system has TO atT > if StopiT) = StopiO) > 0. 

We remark again that the mere presence of non-vanishing TE is not enough to ensure TO. 
Indeed, there are situations for classical states in the sense of [ [27 in which one can encode stable 



classical memory, for instance the toric code in 3D T'^^'^K in which Stop > see also 1351] . In order 



to assess TO, TE must have the full value that it has at zero temperature. In fact, up to half of that 



value can be of classical origin and therefore not related to quantum TO [|25(]. 

We want to discuss the meaning of stable classical and quantum memories from the statistical 
mechanics point of view. The stability of classical memory corresponds to ergodicity breaking. 
The phase space breaks into domains in which the evolution is confined for exponentially long (in 
the system size) times. These domains become disconnected in the thermodynamic limit. One can 
encode classical information in the system by knowing in which domain the system is confined. 
A stable quantum memory means that there is a manifold of metastable quantum states. The 
breaking of ergodicity here is more dramatic, because there is an infinite number of states that are 
disconnected, and that have arbitrarily large overlap between each other. We find that this property 
is accompanied by a particular pattern of long-range entanglement, being it TE or the non triviality 
of entanglement (NT) criterion. We asked ourselves: Why is that? Why should the entanglement 
properties be related in such a strict way to the way the system approaches equilibrium? We think 
that answering this question is crucial to the understanding of the notions of quantum memory, 
ergodicity breaking, and the statistical mechanics of topological phases. 

In this paper, we have proven that the mechanism leading to the destruction of TE is given by 
the confinement-deconfinement transition in the corresponding gauge theory (here the Z2 lattice 
gauge theory) and therefore to a critical temperature Tc for the stability of Stop- Moreover, we 
have proven that the same mechanism is responsible for the transition NT FAC that describes 
topological order in terms of patterns of non trivial long-range entanglement. The confinement- 
deconfinement mechanism though, is also the one responsible for the destruction of quantum (or 
classical) memory [tsij . We can then establish that TO is a property of the way ergodicity is 
broken, or, in information-theoretic terms, is a property of how information encoded in the system 
is resilient. This is important since the entanglement criteria, or other non local order parameters 



|[36ll are hardly experimentally accessible, so we need to find other properties of TO that can 
characterize it. 

This work leaves many open questions, the most natural of all is whether these results extend 
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to general quantum double models and the corresponding discrete gauge theories. Another gen- 
eralization is to investigate this transition in general string-net models 111, [STD, though these do 
not correspond to a gauge theory. Moreover, we have seen why we obtain a sharp phase transi- 
tion when we consider partitions of infinite size. For the finite system, there are exponentially 
suppressed branes that are able to see the topology of the partitions and therefore to give expo- 
nentially suppressed corrections to Stop- With this in mind, we ask ourselves what is the final size 
scaling in order to find the critical behaviour of Stop near (scaling behaviour for the mutual 



information in TO was examined in BSD ). Would the critical exponents for Stop be related to other 
information-theoretic quantities, that we can perhaps measure? Moreover, in the 3D case, we see 
that a TE of classical origin corresponds to the existence of classical memory. Again, we wonder, 
what is the connection? Is there a kind of TE in all the classical phases which host a stable classi- 
cal memory? A fourth question is related to the very interesting class of models whose quantum 
memory at finite temperature is stable for only polynomial (in the system size) lifetimes. These 
models are obtained by coupling with a bosonic system [3^ or with other lon g range interactions 
[|40|] or with Hamiltonians whose excitations have fractal geometric properties [|4ll- l43|l . What kind 
of TO do these models have? What happens to their TE? And what about the Hastings' criterion 
for such systems? 

Along with thermal stability, the question of whether TE in itself is stable in the whole TO 
phase is still an open question in its generality, thought TE is known to be stable in the Z2 lattice 



gauge theory We know that the TO phase at T = is stable because the gap will not close 
for arbitrary perturbations (within a range) [|44]. If TE or NT are the hallmark of TO, one needs to 
prove they are stable too, for the same range. We hope that the results shown in this work can also 
be helpful for this important question. 

To finish, we want to mention that another important form of stability is dynamical stability 
after a quantum quench, that is a sudden dramatic change in the system Hamiltonian. It has been 
shown [|45|] that TE is stable in the quench scenario for some particular quenches in the toric code 
in 2D while quantum memory is not 1461]. Therefore the dynamical scenario is different, and there 
are no results in higher dimensions. Our results in the mapping to the lattice gauge theory may 
prove useful to study also these scenarios. 

All these questions constitute an exciting challenge for the study of novel quantum phases of 
matter, statistical mechanics paradigms, and quantum information theory. 
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